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Abstract 

The author shows the existence of positive solutions to a general 
class of non-autonomous, semilinear elliptic systems in the whole space 
including weighted Hardy-Littlewood-Sobolev type systems. The novel 
method used here implements the classical shooting method enhanced 
by topological degree theory. First, a target map is constructed which 
aims the shooting method, then the non-degeneracy conditions are es- 
tablished which are used to guarantee the continuity of the map needed 
in applying topological degree theory. The existence of zeros for the 
target map guaranteed by degree theory combined with a Liouville type 
theorem for the corresponding Dirichlet problem will show the existence 
of positive solutions to the class of systems. 
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1 Introduction 

A well-known difficulty in studying nonlinear elliptic systems, or any nonlinear sys- 
tem of partial differential equations (PDEs) for that matter, is on developing tools 
useful in their analysis. Such tools are usually limited to specific problems; that is, 
a technique for examining one class of problems may prove ineffective in examining 
other classes of problems. Our main objective of this article is to further develop 
and refine a novel approach — first introduced in [22J — for proving the existence of 
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positive solutions and related properties to higher-order, nonlinear system of ellip- 
tic equations in the whole space. As we shall see below, the class of problems we 
examine will include the well-known Lane-Emden and Hardy-Littlewood-Sobolev 
type systems along with their weighted counterparts as motivating examples. Re- 
markably, the mathematical tools utilized within this framework are more or less 
elementary by themselves. The first main result proves the existence of positive 
solutions, under reasonable assumptions, to the general system 

(-A) ki Ui = fi(\x\,u) iniGl n \{0}, for i = 1,2,..., L. (1.1) 

As will be demonstrated below, proving the existence of positive solutions to this 
system of elliptic PDEs involves reformulating the problem in radial coordinates then 
applying the classical shooting method coupled with a Liouville type theorem for 
the corresponding Dirichlet problem. Specifically, a natural ingredient of the proof 
entails constructing a continuous target map which aims the shooting method. 
Then a fixed point argument via degree theory is invoked to guarantee the existence 
of zeros of this target map which enables us to identify the correct initial shooting 
positions for the shooting method. Combining this with a Liouville type theorem 
will imply the existence of positive, radially symmetric solutions to (jl.ip . 

In establishing our general results, a central and motivating example is the 
generalized weighted system, 

' (-A)iu=^ inIP\{0}, 

' (_A)i« = -j^ inR n \{0}, (L2) 
u, v > in R n . 

Notice that when 7 = 2 and /3i = /?2 = 0, we may consider the whole space and this 
system reduces to the well-known Lane-Emden system, 

f -Au = v g , u > in R n , 

\ -Av = u p , v > in R n , ^ ' 

or more generally the Hardy-Littlewood-Sobolev (HLS) system when 7 > 2 and 
/3i = /3 2 = 0, 

(-A)aii 

{-A)iv 

The Lane-Emden and HLS systems have received much attention in the past few 
decades. For instance, the scalar case was studied in [2JE1IT7], an d similar problems 
have been approached geometrically in [5j[9]. Related systems including its general- 
ized version, the HLS type systems and related problems have been studied as well 



u > in 



v > in 
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0-P2], [H]-[I5], D3H2I]> and [231 [281 EH [35] • When 7 = 2k is an even integer, 
(II. 4p is equivalent to the integral system 

f v q (y) 

u(x) = / A^dy, u > in M n , 

Jm n |x — y\ n t 

u(z) = / , UP{V l_ dy, OO in R n , ^ 

Jr« \ x — y\ ' 

< p, q < 00, < 7 < n, 

in the sense that a solution of one system, multiplied by a suitable constant if 
necessary, is also a solution of the other when p, q > 1. Hence, systems (|1.4j) and 
(jl.5p are both referred to as the HLS system. Now when studying the HLS system, 
the exponents p, q, and 7 play an essential role in its behavior. More precisely, there 
are three important cases to be considered: for the positive constants < p, q < 00 
and 7 = 2k (k G N), 

1 1 n-2k 

w — p + TTq' > ~ — ' 

1 1 n-2k 
lb) 1- = , and 

v ; 1+p 1 + q n 

1 1 n-2k 

w 1 +p 1+q n 

The HLS system is said to subcritical if (a), critical if (b), and supercritical 

if (c). In the case of (jl.3p . the famous Lane-Emden conjecture — an analogue to 
the celebrated result of Gidas and Spruck [T7] in the scalar case — states that this 
elliptic system with subcritical exponents has no classical solution. This has been 
completely settled for radial solutions [271 US] , for dimensions n < 4 [301 EH EZ] , and 
for n > 4 under certain subregions of subcritical exponents [31 [161 EH E21 EH EE]- 
With the help of the moving planes method in integral forms, the work in [13j— 
when combined with the non-existence results in [26] — provides a partial resolution 
of this conjecture as well. 



In the critical case, the integral system (jl.5p is the Euler-Lagrange equations 
of the fundamental Hardy-Littlewood-Sobolev inequality. In this case, we call a 
solution of the above system ground state if it 'optimizes' the HLS inequality. 
Recall that the HLS inequality (see [El [25J [36] ) states that 

f . {x)9{ y} dxdy<C sXn \\f\\ r \\g\\ s (1.6) 
\x — y\ A 

where 0<A<n, l<s,r<oo, I + i + ^ = 2, / G L r (M n ), and g G L s (W l ). To 
find the best constant in the HLS inequality, one maximizes the HLS functional: 

J(f.9)-[ I f -^d*iy (1.7) 
Jm n Jm n f — y\ 
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under the constraints 



\fJ\ 



1. Let p = q = ~r and with a suitable 



scaling such as u = c\f r 1 and v 



C29 



s-l 



the Euler-Lagrange equations are 



precisely the system of integral equations in (jl.5p . Here, u £ L' p+1 and i> S L 9+1 are 
equivalent to / € L r and g € L s and the exponents, < p, g < 00, satisfy 



1 



+ 



1 



n 



7 



p + 1 q + 1 



Moreover, Lieb proved in |25j the existence of positive solutions to (jl.5p which 
maximize the corresponding functionals J{f,g) in the class of it € L p+1 and u € 
In other words, there exist maximizers to the Euler-Lagrange equations under 
critical exponents, therefore showing the existence of ground state solutions for the 
HLS system. In addition, Hardy and Littlewood also introduced the following double 
weighted inequality which was later generalized by Stein and Weiss [39] : 



f(x)g(y) 



— y\"\y 



dxdy < C a ,p iS ,\,n 



\9\\ 



(1.8) 



where a + j3 > 0, 



1 X a 
1 < - < 1 

r n n 



1 ,11 A + a + /3 
-, and — I 1 

r r s n 



To find the best constant in the double weighted inequality, one maximizes the 
associated functional: 



J(f,g) 



f(x)g(y) 



\xr\x 



y\ x \y\ 



rdxdy. 



As was shown above, one can show that the corresponding Euler-Lagrange equations 
are the system of integral equations: 



1 



x 



v q (y) 



l\P\x - y 
uP(y) 

y\ x 



jdy, u > in 
dy, v > in '. 



(1.9) 



u(x) 
v(x) 

pt J«. n \v\~\ x - y\" 

where 0<p,q<oo, 0<A<n, a < J^. < an d ^ + ^ = Interest- 
ingly, the authors in [llj considered a weighted HLS system. Namely, they proved, 
under appropriate mild conditions, the uniqueness of solutions to the singular non- 
linear system, 



-A(\x\ a u) 
-A(\x\Pv) 



m 



n \{0}, 
in M n \{0}, 



11.10) 



and classified all the solutions for the case a = (3 and p = q, thus obtaining the best 
constant in the corresponding weighted HLS inequality. 
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Let us remark on the case of supercritical exponents for the HLS system in 
relation to the work in this article. As a simple illustration, let u = v and p = q in 
(|1.5p to obtain the following scalar integral equation, 

r u p (v) 

u{x) = / {J T^zdy, 00 in R n , (1.11) 

jR n \ x — y\ 1 

with the corresponding partial differential equation 

(-A) k u(x) = u p (x), 2k < n, u > in R n . (1.12) 



As before, both the integral and differential equations are called supercritical if 
P > ^r^i critical if p = and subcritical if p < In the supercritical case 

with k = 1, the shooting method can be successfully applied to (|1.12p . however, 
much difficulty arises even in the scalar case for k > 2. In the results of this pa- 
per, we circumvent these difficulties by enhancing the shooting method via degree 
theory in an interesty way and the approach used applies to systems as well. So 
the primary objective of this paper is to further develop our framework to handle 
even more general systems such as the weighted system (II. 2\i especially since such 
existence results are not so well developed for these problems. Hence, we shall deter- 
mine the conditions on such systems which allow us to prove existence of solutions 
using our technique, and in doing so, we demonstrate how to handle even the case 
of (II. 2p which is not included in the results of [22\. In addition to the Liouville type 
results, the difficulty in implementing our technique lies in determining the sufficient 
conditions which guarantee the continuity of the target map. This difficulty moti- 
vates our consideration of non-degeneracy conditions. Specifically, we introduce 
non-degeneracy conditions each geared to handle systems such as (|1.2|) with varying 
exponents and weights. 

The rest of this manuscript is structured as follows. In section [2j we introduce 
some preliminary definitions and the precise statements of our main results. Section 
[3] gives the proofs of the first existence theorem concerning the general system . 
In section U we prove the existence theorems dealing with the weighted system 
(II. 2p . In order to prove this theorem, a Liouville type theorem for this system is 
required; thus such a non-existence result for the corresponding Dirichlet problem 
is also provided in this section. 

2 Preliminaries and Main Results 

Consider the system 

-A) k *u i = f i (\x\,u), inM«\{0}, 

Ui > 0, in R n , for i = 1,2,... ,L. 1 ' ' 

with the following assumptions. From this point on, we will always assume ki > 1 
and F(\x\,u) = (/i(|x|, n), f2(\x\,u), . . . , /l(|x|, u)) satisfies the following conditions: 
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is a continuous vector- valued map, 

(b) F(\x\,u) is locally Lipschitz continuous in the second argument uniformly in the 
interior of M + x . 

Non-degeneracy condition 1. Given a v G dW}, let I® and l£ contain the 
components j G {1,2,..., L} of v such that = and v j > 0, respectively. Then 
there are constants A = X(v) > and /3 > —2 and a 5 = 5(v) such that if \v — w\ < 5 
then 



iexg 

Remark 2.1. In £/ie case where k{ = 1 /or all 1 < i < L, we shall always assume 
that fi (\x\,u) = fi (u) for some 1 < io < L. 

Non-degeneracy condition 2. Suppose that ki = 1 for all i^(|a;|,u) = F(v), 
and ^(v) 7^ whenever v > in (|2.ip . Given any two real numbers < m < M, 
we have for v = (yi,V2, ... , vl) € such that 

Vi k < for = 1, 2, . . . , j and m < Vi k < M for A; = j + 1, . . . , L, 

where . . . , il} is any permutation of {1,2, ...,£}, there exists a constant 

C m ,M > such that 



Moreover, if the right-hand side of (|2.2ft is zero, then the left-hand side must also 
vanish simultaneously. 

Systems of the form (|1.1|) satisfying conditions (a)-(b) and the non-degeneracy 
'condition 1' (or 'condition 2') are said to be non-degenerate 'type I' (or 'type 
II'). The weighted HLS system is an example of a non-degenerate type I system and 
the system, 



is an example of a non-degenerate type II system provided p > t > 0, and q > s > 0, 
but also notice that the latter is not of type I. 

The first theorem presented in this paper illustrates how the existence of solu- 
tions for non-degenerate systems (jl.ip follows from the non-existence of solutions to 
the corresponding Dirichlet problem, 




max 

j+l<k<L 




(2.2) 



k=l 




(2.3) 




: 'Ui = fi{\x\,u) 
Ui > 



in £ fl (0)\{0}, 
in B R (0), 
Ui = on dB R (0), i = 1,2,... 



< 



(2.4) 



k Ui = -Am 



... = (-A) k *- 1 



L 
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for all R > 0. Here -Br(O) denotes the open ball of radius R centered at the origin 
and sometimes we use B R instead for brevity. We will be mainly concerned with 
type I systems in this paper, but the authors in |24j obtained analogous results for 
unweighted, non-degenerate systems of type II. Namely, as a motivating example, 
one of the results in that work proved the following: 

Theorem. The non- degenerate system (|2.3p has a solution of class C 2 (IR n ) provided 
that q >t >0, p > s > 0, and 



1 



+ 



1 



< 



n 



l + q ' l+p ~ 2 
Now we are ready to state our main results. 



Theorem 1. The non- degenerate type I system (|2.ip admits a radially symmetric 
classical solution provided that (12. 4ft admits no radially symmetric classical solution 
for all R>0. 



This theorem has the following consequence. 
Theorem 2. The system 

(-A) k u 
(-A) k v 



x 



\S 2 



in R n \{0}, 
in W l \{0}, 



u, v > 



(2.5) 



admits a solution of class C 2k (M. n \{0}) provided that 

n - P\ n — P2 



l+q + l+p 

where < p, q < 00 and /3i, /?2 € ( — 00, 2). 



<n-2k, 



In addition, we have the following Liouville type result which will be important 
in proving Theorem [2j Notice that we can show non-existence of solutions to an 
even more general system than the weighted HLS system. 



Theorem 3. The system, 



-A) k u 



(-A) k v 



U V 



ytyj) 



m 5 fl (0)\{0}, 
in B R {Q)\{0}, 



u 

V 



u, v > in -Br(O), 

-Au = ■■■ = (-A) k ~ l u = ondB R (0), 



(2.6) 



-Av 



(-A) 



k-l„ 



ondB R (0), 
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admits no radially symmetric solution of class C 2k (Bf((0)\{0}) n C 2k 1 (Br(0)) for 
any R > provided that 

< n — Ik 2.7 

l+q 1+p ~ y ' 

where s,t,p,q > and Pi, 02 £ ( — oo, n )- 



3 Proof of Theorem U 



In order to prove the first two theorems we must introduce several key ideas and 
lemmas. As mentioned earlier, the proof centers on constructing a map which aims 
the shooting method. A crucial step will be to show the continuity of this map, and 
the non-degeneracy conditions are exactly what is needed in guaranteeing continuity. 



For i = 1, 2, . . . , L set iOj j = (—A) 



3-1 



U;. 



1 < 3 < k( and consider the system 



( - AlVi t l = w it2 , 

- &Wi,2 = Wi,3, 

- Aw itkz = fi(\x\,wi t i, io 2) i, • • • ,wl,i) in M n \{0}, 
Wi,l,Wi,2, ■ ■ -,Wi,ki > in R n , 

where i = 1, 2, . . . , L. 



(3.1) 



Solutions to (13. lj) are clearly solutions to (12.1j) , so it will suffice to show the existence 
of solutions to f)3. 1 1) instead. We can express the above system into the form: 



f - Aioi = fi(r, w), -Aw 2 = f2{r,w), 

- Aw 3 = f A {r,w),. . . ,-Aw L -i = f L -i(r,w), 

- Aw L = f L (r,w) in 
k 101,102, • • • ,wl > in 



r\{o}, 



(3.2) 



where we still use L to represent a generic positive integer. We shall work with (|3.2p 
instead when proving Theorem HJ but note that the non-degeneracy condition 1 still 
holds true for this new system. Now let us define the aforementioned target map. 
For any strictly positive initial value a = (ai, «2, • • • , cth), consider the IVP 



u>i ( r ) + w i( r ) 



-fi(r,w(r)), 



[ w-(0) = 0, i0j(0) = ai for i = 1,2, 
Clearly (|3.3p is equivalent to (|3.2p in radial coordinates. 



L. 



(3.3) 
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Definition 3.1. Define the target map ip : h-> as follows. For a G mi(R^), 
i/ie interior o/IR+, 

f'aj VK ) = ^(^0) where ro is i/ie smallest such r for which Wi (r) = for some 
l<i <L, 

(b) otherwise, if no such ro exists, then ip{a) = lirn r _ >tX) u;(r). 

(c) Moreover, tp = Identity on cM^ . 

Remark 3.2. We may think of (a) as the case when the solution hits the wall for 
the first time and (b) is the case where it never hits the wall. Observe also that ip is 
equivalent to the identity map on the wall. This property is crucial when we apply 
the tools from topological degree theory. 

The next lemma is a standard result from Brouwer topological degree theory 
and can be found in various literature, see [1] and [29] for instance. 

Lemma 3.3 (Dependence on boundary values). Let U C M n be a bounded open 
set and f,g : U i-> W 1 are continuous maps. Suppose that f = g on dU and 
a ^ f(dU) = g{dU), then degree(f,U,a) = degree(g,U,a). 

One may recall the important property that if degree(f, U, a) 7^ 0, then there exists 
a point x G U such that f(x) = a. 

Lemma 3.4. The target map ip : 1— > cM;^ is continuous. 

Proof of Lemma 1. Choose any a G and fix an e > 0. Without loss of generality, 
there are two cases to consider: 

(1) a G dR^, 

(2) a G int(R+) and ip(a) = w(ro,a), 

We only consider the two cases since the remaining case where a G int(W+) with 
w(r, a) > for all r > says that we have found an initial shooting position a that 
gives us a desired global solution as stated in the main theorem. 

Case (1): By definition we have that if) (at) = a and fix e > 0. We can find a 5i > 
such that for \a — a\ < 5i = 5\(a), the non-degeneracy condition 1 implies that 

^2 fj(r,a) > \(a)rP for all r < 1. 
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Then we can find 62 > such that \a — w(r,a)\ < 8% for r < 82 and \a — a\ < 82- 
Let 

W°(r,a) := Wj(r,a). 
jell 



Then the non-degeneracy condition 1 and (|3.3h imply that 

dr \ dr J 
Integrating this twice with respect to r yields 



u'eJS- 



From this we can find 5 > sufficiently small with <5 < e such that W°(r a , a) = 
for some r a < 62 and 

|?/>(a) — ^(a)| = \ip(oc) — a\ < \w(r a , a) — a\ < 5 < e 

whenever \a — a\ < 5. 

Case (2): Since the source terms fi are non-negative, u' io (ro,a) < by a direct 
computation or simply by Hopf's Lemma. This transversality condition along with 
the ODE stability imply that for a sufficiently close to a, the solution to this per- 
turbed IVP must hit the wall and ip(a) must be close to ip(a). 

This completes the proof that tp is continuous at a € M.+ . □ 
Lemma 3.5. For every a > 0, there exists an a a € A a where 

A a := |q E I Yl 

such that ip(a a ) = 0. 

Proof of Lemma 2. Define the set B a as follows 



«i = a 



B a := < a E dRi 



^ ctj < a > ■ 
i=i J 



It follows that ip maps A a into B a due to the non-increasing property of solutions. 
Now define the continuous map (j) : B a — > A a by 



(a) = a + 
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with continuous inverse 



: A n 



B n denned 



<fT\a) 



a 



mm on 



Set rj = 4> o ip : A a — > A a . Then r/ is continuous on A a and is equivalent to the 
identity map on the boundary of A a . By Lemma |3.3| the index of the map satisfies 
degree(r], A a ,a) = degree{Identity, A a ,a) = 1 7^ for any interior point of A a . So 
r] is onto, and thus ip is onto. Then there exists an a a £ A a such that ip{a a ) = 0. □ 



Proof of Theorem [p Let w = w(r) be a solution of f|3.3[) with initial position 
w(0) = a a as guaranteed by Lemma [3. 5 1 for fixed a. This solution must never hit the 
wall. If this was the case, then there would be a smallest finite value r = tq such that 
w(fo) = ^{oia) = 0. But this would imply that w = w{\x\) is a radially symmetric 
solution to (|2.4p with R = ro, which is a contradiction. Hence, w = w (|x|) must be 
a radially symmetric solution to (|2.ip . 



Remark 3.6. In the proof of Theorem {J\ we are implicitly using the fact that the 
Dirichlet problem (12. 4h is equivalent, under classical solutions, to the Dirichlet prob- 
lem, 



- Awi,i = w i>2 , 

- Aw ji4r i = w kt , 

- Aw itki = fi(\x\, Wl > l,W2,l, 

Wi,2, ■ ■ ■ ,Wi,ki > 
^ W^i = Wi,2, ■ ■ ■ ,m,ki = 



,w L ,x) in 5 fl (0)\{0}, 
in B R {0), 
on dB R {0), 



where Wij = (— A) J_1 Uj for j = 1, 2, . . . , ki and i = 1,2, ... ,L. This equivalence 
follows from an inductive argument by exploiting the non-negative property of the 
source term F, the boundary conditions, and the strong maximum principle. We 
shall revisit this property in the proof of Theorem 
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4 The Weighted System 



The system (|1.2p can be reduced into the system 



Awi = 
Aw k = 
Aw k+ i 
Aw 2 k = 



101,102, 



W2, ■ ■ 

W l W k+l 

\x\h ' 

= W k+ 2, ■ ■ 
. W k+l W l 

,W 2 k>0 



-Aw k - 



Wk, 



-Aw 2 k- 



W2k, 



(4.1) 



m 



m 



l \{0}, 



where W\ := u and Wk+i '■= v. Let Br := Br(0). From Theorem[T]it suffices to show 
the non-existence of solutions to the system (14, 2D given below but under Dirichlet 
boundary conditions on balls. That is, we will show that the system 

- Awi = U>2, 



Awk-i 
Awk+i 



w k ,-Aw k 

Wk+2, 



(4.2) 



AlV2k-l = W 2 k, -Aw 2 k 



t P 
W k+l W \ 



X 



$1 



in B R \{0} 



wi,w 2 ,...,w 2 k > 

\Wl=W 2 = ■■■ = W 2 k 







111 

on 



Br 
dB R 



admits no solution of class C 2 (Br\{0}) n C 1 (Br) for any R > 0. Observe that the 
usual Liouville type results such as those found in [27] and [31] do not apply since the 
source terms in (14. 2p do not follow from a potential function so that a variational 
approach cannot be used. Therefore a result for the non-existence concerning this 
weighted system, Theorem [3l is provided in order to bypass this issue. First, we 
need the following lemma. 

Lemma 4.1. Let Br again denote the ball in W 1 of radius R centered at the origin, 
and let wj (J = 1,2, ... , 2k) solve |^.j?| ). Then 



u s v g+l 
~\xW 



dx 



yt u P+l 



B R \ x 

/ Vwj ■ Vw 2 k+i-j dx 

J B R 



dx 



12 



/ Wj+iW2k+i-j dx =: E. (4.3) 



I Br 

Proof. To prove this lemma, multiply the 2/c-th equation in (j4.2j) by w\ then in- 
tegrate over Br. The repeated application of integration by parts along with the 
boundary conditions yield 



t „,j>+i 

dx = I \7w\ ■ Vu>2fc dx 



w k+i w i 



B R N & J Br 



W 2 k^Wl dx = I W 2 kW2 
Br J Br 

W2^W2k-l dx = Vu>2 ■ Vu^fc-l dx 

Br J Br 

w 2 k-iAw 2 dx = / w 2 k-iw?>dx 

Br JBr 



Vw k ■ Vw k+ i dx = - w k+1 Aw k dx 

Br J Br 

s 9+1 
W l W k+l , 

— i — rs — dx. 

Br 



□ 



Remark 4.2. Let us be more precise in the calculations found in our proof of Lemma 
4-l\ since we will be using similar calculations below. For instance, when we multi- 
ply, say, the 2k-th equation —Aw2k = Wfc+i^f with w\ then integrate over the ball 
Br, this should be understood implicitly in the following way. We integrate over 
Bn\B t (0) for < e < R and use an integration by parts to obtain 



t „„p+i 

dx = — I w\Aw2kdx 



u n -u,^\ vA 132 ' J«\ "<(()) 

w\ ds + [ Vwi -Vw2kdx, 

8B t {0) OV Jb r \B,(0) 

where v is the inward unit normal vector along dB e (0). Taking the limit as e tends 
to zero, the surface integral will vanish since the Wi's are of the class C 1 (Br(0)) 
and are therefore uniformly bounded on Br along with their first- order derivatives. 
Then we obtain 

w p+1 f 
t 1 ,. 1 — dx = Vwx ■ Vu?2fc dx. 

'Br M & JBr 

All such calculations including those found in the proof of Theorem \3\ below should 
be understood in this way. 
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Proof of Theorem [3l The equivalence of (|2.6p and f)4.2[) under classical solutions 
follows from the maximum principle. To see this, first observe that if u = w\ and 
t> = Wk+i where w^s satisfy (|4.2p . then -u and t> must also satisfy (|2.6p . Now suppose 
u and v satisfy (12. 6p and let = (— A) l_1 u for i = 1, 2, . . . , k and = (— A) l ~ k ~ 1 v 
for i = A; + 1, k + 2, . . . , 2k. We only need to show the super polyharmonic property: 
Wj's > in Br. Since it is already given that w\ and Wk+i are positive in Br, we 
have that —Aw^, —Aw2k > in Br from (|4.2p . The Dirichlet boundary conditions 
along with the strong maximum principle imply that wt,W2k > in Br. Then it 
follows from (|4.2p that -Akj^, - Atti 2 i_i > in Br. We can repeat the same 
argument to show that Wk-\,W2k-i > in Br. In fact, we may inductively repeat 
this argument to show the remaining w^s are positive in Br, thus completing our 
verification that the two systems are equivalent. So with this in mind, it will suffice 
to show that system (I4.2p admits no solution of class C 2 (Br\{0}) n C 1 {Br) under 
the constraint (|2.7p with s,t,p,q > 0. 

For j = 2,3, ... ,k — 1, multiply the j-th equation in (|4.2p by x ■ S/w2k+i-j, 
integrate over Br, then integrate by parts to obtain 

-f ^ ^2 fc +W {x . n)ds+ f Vwj . Vw 2k+1 ^ dx (4.4) 
JdB R an on j Br 

+ / x ■ w jtXi V(w2k+x-j)xi dx = w j+ i(x ■ Vw 2 k+i-j) dx, 

JBr JBr 

where n is the outward pointing unit normal vector. Multiply the (2k + 1 — j)— th 
equation in (|4.2p by x-Vwj and integrate over Br and perform analogous calculations 
as was done in obtaining (|4.4p . Then summing the resulting equation with (|4.4|) and 
using the identity, 

X ■ Wj tXi V(w2k+l-j)xi + X ■ W 2 k+l-j,x i 'V(Wj) Xi dx 

Br 



/ x ■ V(Vio,- • Vw 2 k+i-j) dx 

JB n 



dwj dw 2k +i - 
b r " ~ J ' ~' K ^ l ~' J ~" ' J aBR dn dn 



n j Vwj • Vu>2fc+i-j dx + j — ^q^ 1 ° (x ■ n) ds, 



we obtain 



dw 1 dw 2 k + i- ] {x n) dg + (2 _ n) j y W2k dx (4.5) 

9Br dn dn J Br 

w j+1 (x ■ Vw 2k+ i-j) + w 2 k+2-j(x ■ Vwj) dx. 

Br 



Now multiply the 2/j-th equation in (|4.2p by x ■ Vwi and integrate over Br to obtain 
(x ■ Vu>i)Au> 2 fc dx = \ (x ■ Vwi) k+ \ 1 dx . 

Br " JBr \X\PZ 
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Let us calculate I\ and 12- Using integration by parts, 

h = - [ ^p-^^{x ■ n) ds + [ Vw 1 -Vw 2 kdx+ [ x^^- ( f Wl ) <l.r. 
JdB R On dn J Br J Br dxj \dxjdxi 



and 



n-fh [ w Ui< +1 dx t [ Oi +1 fx Vm l)dx 
Now multiply the first equation by x ■ Vw2k and integrate over B B to obtain 



/ (x ■ Vui2fc)Auii dx = (x ■ Vw 2 k) 

JB R JBr, 



W2 dx . 

Br 



:=JJl :=Ih 

We use integration by parts to rewrite II\ as follows. 

IIx = -[ 9 ™ 2k (x ■n)ds+ [ Vw! ■ Vw 2 k dx + / Xi^p- ( f ^ 2fc ) dx. 
JdB R on dn J Br J Br dxj \dxjdxij 

By summing together the two equations I\ = I2 and II\ = II2 and using the fact 
that 

/ x ■ V(Vwi ■ Vw2k) dx = [ ^P-^p-(x ■ n) ds — n [ Vwi • Vit^fc dx, 
Jb r ' JdB R on dn J Br 

we obtain the identity 

(2 — n) f Vw2k ■ Vwi dx + — — — / k ~!~ 1 l a 1 — dx 
Jb r ' 1+P Jb r fp 

= / ^p^-Tp-ix ■ n) ds + / w 2 (x ■ Vw 2 k)dx (4.6) 
JdB R on dn J Br 

h w k+i 



1 +PJb r 



\.r 



1/32 



(x ■ S/wk+i) dx. 



Multiply the fe-th and (A: + l)-th equations in (|4.2p by x ■ Vwk+i and x ■ S/iVk, 
respectively, and integrate over Br. Using similar calculations to those used in 
deriving (|4.6p . we obtain 



(2 - n) / \/w k ■ vwk+i dx + —— / „ dx 

Jb r 1 + 9 Jb r FP 

dw k dw k+ i 



9Br dn dn J Br 
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x-n)ds+ I Wk+2(x ■ Vwifc) dx (4.7) 



9+l„„s-l 

[x ■ Vw\) dx. 



w k+i w i 



Observe also that integrating by parts and using the boundary conditions, we see 
that 

/ x ■ (wj + {Vw 2 k+l-j + u>2fc+i-jVu>j+i) dx = x- V(w j+ iw 2 k+i~j) dx 
Jb r Jb r 

= - n I Wj + iw 2 k+i-j dx. 



Using this identity and summing (|4.5p over j = 2, 3, . . . , k — 1 along with (|4.6p and 
(|4.7p . we see that 



(2 — n) \^ / Vwj • Vt«9/s+i-j dx 



1 +Q J Br \ X \ Pl ^ 1 +P J Br \X\ 02 



n / Wj+iW2k+i-j dx 
j=l J Br 

f^JdBRdn dn [X n)dS l + q J BR [X VWl)dX 



1 +PJb r \x 



' x ■ Vwk+i) dx. 



Observe that the right hand side of this inequality must be strictly positive by the 
non-increasing property of the positive radial solutions. Hence, Lemma 14.11 implies 
that 

n — /3i n - /3 2 
+ q 1+p 

In other words, we have 



( , , n — pi n — Ho , . ~i 

<^ k(2 - n + - — — + - — — + (k - l)n \ ■ E > 0. 

I 1 + a 1 + p ) 



1+q l+p 

but this contradicts with (12.71). □ 



Proof of Theorem [2l The continuity of the target map guaranteed by the non- 
degeneracy condition along with the Liouville type result of Theorem [3] with s = 
t = will imply the desired result as a consequence of Theorem [TJ □ 
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Remark 4.3. In the autonomous case where F(\x\,u) = F(u) with the additional 
assumption that u = whenever F(u) = ; the definition of the target map implies 
that the positive solutions will vanish at infinity. Namely, basic elliptic theory implies 
that F(if)(a a )) = 0, so that ip(a a ) = from this extra condition. In other words, the 
global positive solutions given by our results must exhibit the following asymptotic 
behavior: 

Ui — > uniformly as \x\ — > oo, for i = 1, 2, . . . , L. (4.8) 

The desire for similar types of asymptotic behavior is precisely why certain assump- 
tions where placed on F in the non- degeneracy conditions; that is, we required such 
conditions so that in the case where the positive radial solution never hits the wall, 
at least one component of the solution decays uniformly to zero and the target map 
takes values on the wall. This asymptotic property was key in showing the target 
map was onto, as well. 
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